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ABSTRACT 

We consider the thermal properties of cold, dense clouds of molecular hydrogen and atomic helium. For cloud 
masses below 10~ 1,7 Mq, the internal pressure is sufficient to permit the existence of particles of solid or liq- 
uid hydrogen at temperatures above the microwave background temperature. Optically thin thermal continuum 
emission by these particles can balance cosmic-ray heating of the cloud, leading to equilibria which are thermally 
stable even though the heating rate is independent of cloud temperature. For the Galaxy, the known heating rate 
in the disk sets a minimum mass of order 10 _6 Mq necessary for survival. Clouds of this type may in principle 
comprise most of the dark matter in the Galactic halo. However, we caution that the equilibria do not exist at 
redshifts z >1 when the temperature of the microwave background was substantially larger than its current value; 
the formation and survival of such clouds to the present epoch therefore remain open questions. 

Subject headings: ISM: clouds — galaxies: halos — dark matter 



1. INTRODUCTION 

Walker & Wardle (1998) showed that a population of neu- 
tral, AU-sized clouds in the Galactic halo could be responsible 
for the "Extreme Scattering Events" (ESEs) observed in the ra- 
dio flux towards several quasars (Fiedler et al. 1987, 1994). 
In this model the cloud surfaces are exposed to UV radiation 
from hot stars in the Galactic disk, producing a photo-ionised 
wind. When one of these clouds crosses the line of sight to 
a compact radio source, the flux varies as a result of refrac- 
tion by the ionised gas (cf. Henriksen & Widrow 1995). This 
model explains the observed flux variations quite naturally; but 
if the clouds are self-gravitating, then the ESE event rate im- 
plies that the cloud population comprises a significant fraction 
' of the Galaxy's mass. 

This halo cloud population cannot contain much dust mixed 
with the gas as this would lead to optical extinction events of 
distant stars: either the clouds have extremely low metallicity, 
• or any dust grains have sedimented to the cloud centre. Given 
this, several factors make the clouds difficult to detect (Pfen- 
niger, Combes & Martinet 1994): cold molecular hydrogen is, 
by and large, invisible; the clouds are small; they are transpar- 
ent in most regions of the electromagnetic spectrum; and they 
cover a small fraction of the sky. The clouds are not sufficiently 
compact to cause gravitational lensing towards the LMC, al- 
though Draine (1998) has shown that there is substantial optical 
refraction by the neutral gas, so that microlensing experiments 
(Paczyhski 1996) already place useful constraints on the prop- 
erties of low-mass halo clouds. 

Given that this hypothesised cloud population does not vio- 
late observational constraints, the primary issues that need to 
be addressed are theoretical: (i) how and when did these clouds 
form? and (ii) how do they resist gravitational collapse? The 
second of these is addressed in this Letter. 

We begin by writing down equations describing a simple 
"one-zone" model of a cloud, characterised by a single temper- 
ature and pressure (fQ), and show that particles of solid H2 may 
exist in the clouds (cf. Pfenniger & Combes 1994). At tem- 
peratures above the microwave background temperature, these 
particles cool the cloud by thermal continuum radiation, admit- 



ting equilibria in which this cooling balances heating by cos- 
mic rays. In §[}J we show that (for optically thin emission) these 
equilibria are thermally stable: if the cloud contracts the coolant 
is destroyed by the increase in temperature, and the power de- 
posited by cosmic rays causes the cloud to expand and the tem- 
perature to return to its original value. We conclude that, within 
the context of our one-zone model, the viable mass range for 
Galactic clouds is 10~ 6 -10~ 17 Mq. 

2. CLOUD MODEL 

Virial equilibrium implies that for a self-gravitating cloud 
characterised by mass M, temperature T, and radius R, 



R w GMfi/kT 



(1) 



where /i is the mean molecular weight. The pressure in the 
cloud can be related to the temperature upon noting that P ~ 
GM 2 /R A , yielding 



P = 



G 3 AP 



kTY 



(2) 



where q depends on the cloud's structure. For poly tropes q rises 
monotonically from 9.2 to 40 as the polytropic index runs from 
3/2 to 9/2, so we adopt q = 20. 

At sufficiently high pressures a fraction x of the molecular 
hydrogen assumes solid (or liquid) form. Then fi = (1 — x + 
2y)m/ (1—x + y), where m is the mass of an H 2 molecule, and 
y w 1/6 is the abundance ratio He:H2 by number. Neglecting 
the temperature difference between the phases, in equilibrium 
the partial pressure of H2 equals the saturated vapour pressure, 
i.e. 
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(valid for < x < 1), where kT v is the heat of vapourisation 
for H 2 (Phinney 1985). With T v = 91.5 K, the vapour pressure 
given by the RHS of eq. (g) is within 20 % of the available 
experimental data (Souers 1986). 

Hydrogen grains can cool the gas in a manner similar to dust 
grains in molecular clouds: the gas cools via collisions with 
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FIG. 1 . — Top: cooling by rotational transitions of H2, HD and LiH in the gas phase, and thermal continuum emission from particles of solid H2 for a 10 — 3 Mq 
cloud in virial equilibrium. Bottom: solid H2 cooling rates for cloud masses 10~ 3 , 1(T 5 and 1O- 7 M . Dotted portion of the IO^Mq curve indicates where 
the coolant is liquid H2. The dashed line in both plots indicates the cooling rate that would give a Kelvin-Helmholtz time scale of 10 Gyr. 
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slightly colder solid particles, which in turn cool by thermal 
continuum emission. To calculate the cooling by solid H 2 , first 
consider the net power radiated by a single particle (we employ 
an 'escape probability' formulation of radiative transfer): 



4cr 



/C(T)T 4 
V 1 + t 



C(T b )T£ 

l + n 



(4) 



where C (T) is the Planck-mean absorption cross-section, T\, is 
the cosmic microwave background temperature, and r and r& = 
tC(Tj,)/C(T) are the Planck-mean optical depth of the cloud 
to thermal radiation characterized by T and T b respectively. As- 
suming that the particle size is <C A (~ 0.1 cm at the temper- 
atures of relevance here), we may write C(T) = C m (T)m s , 
m s being the particle mass, and for spherical grains we have 
(Draine & Lee 1984) 
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where p s — 0.087 gcm~ 3 is the density of solid H 2 (Souers 
1986), the complex dielectric function of the solid is e-y + ze 2 , 
and we have assumed that e 2 = A 2 / A as expected at low fre- 
quencies. The net cooling rate per unit mass of cloud material 
(gas and solid) is then 
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where n = r (T b /T) 2 , and 
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To evaluate A, we require optical constants for solid H 2 in 
the microwave. The particles are expected to be almost pure 
para-hydrogen as an ortho-para mixture of the solid relaxes to 
para (J = 0) form in a few days (Souers 1986). The low 
frequency value of s% for para-hydrogen has been measured 
(Souers 1986) as E\ ~ 1.25; the low frequency limit of £ 2 is 
less certain. Jochemsen et al. (1978) measured the extinction 
coefficient of a single crystal of solid para-hydrogen in the re- 
gion of interest (A ~ 0.1 cm), but could not determine whether 
this continuum extinction was due to absorption or scattering 
within the crystal. Because these measurements do not con- 
form to the anticipated low-frequency behaviour (oc 1/A), and 
absorption bands are not expected below the S(0) line, it is 
likely that the absorption of pure crystalline H 2 is much smaller 
than the measured extinction, and we can only infer a limit: 
e 2 <^ 1.8 x 10~ 3 . However, the low-frequency absorption of 
solid H 2 grains could be strongly enhanced by impurity species 
and lattice defects. For the purposes of this paper we adopt 
£2 = A 2 /A and A 2 = 10~ 4 cm. Within the confines of the 
model, this assumption represents one of our main areas of un- 
certainty. 

For a given cloud mass, the fraction of H 2 in the solid phase 
can be determined from T using eqs (Q) and This allows 
the cooling by solid particles to be calculated as a function 
of T. The upper panel of Fig [j] illustrates this for a cloud 
of mass 10 _3 Mq. For comparison, the cooling contributed 
by rotational lines of gas-phase H 2 , HD and LiH is also plot- 
ted; energies and A-values from Turner, Kirby-Docken & Dal- 
garno (1977), Abgrall, Roueff & Viala (1982), and Gianturco 



et al. (1996). The adopted deuterium and LiH abundances are 
3 x 10~ 5 , and 1.2 x 10~ 10 respectively (Schramm & Turner 
1998). We employ an escape-probability formulation of radia- 
tive transfer, in which the optically-thin cooling rates are di- 
vided by (1 + t) where r is the optical depth at line centre (ro- 
tational transitions) or the Planck-mean (continuum emission). 

There is a critical temperature T c at which H 2 in the cloud 
lies on the border between the solid and gaseous phases, i.e. 
the partial pressure of H 2 is equal to its saturated vapour pres- 
sure. From this point x and A increase precipitously as T is re- 
duced, until t ~ l and the cooling is then roughly black-body. 
A then increases more slowly, peaking when x ~ 1/2 and sub- 
sequently dropping to zero as the temperature approaches that 
of the microwave background. At lower temperatures the cloud 
is heated by the background radiation. 

The solid H 2 cooling curves for cloud masses of 10~ 3 , 10~ 5 
and 10~ 7 Mq are compared in the lower panel of figure [T] 
Decreasing the cloud mass has two consequences: the criti- 
cal temperature T c decreases, whereas the maximum value of 
A decreases (because in this circumstance the emission is op- 
tically thick, and cloud surface area is proportional to M 2 at 
a given temperature). For cloud masses < 10 _7 Mq, T c is 
above the H 2 triple point (13.8 K) and liquid droplets of H 2 
form instead. This does not qualitatively affect the calculations 
as the density and saturated vapour pressure of the liquid are 
within 50% of those of the solid for T £ 20 K (Souers 1986), 
and the optical properties are similar (in the sense that e 2 is 
small and uncertain), thus in Fig. [Fj we continue the cooling 
curve for a 10~ 7 Mq cloud above the triple point as a dot- 
ted curve. Thermally stable solutions do not exist for masses 
<^ 10 _7 ' 5 Mq, as the partial pressure of H 2 exceeds the satu- 
rated vapour pressure unless x > 0.5 (see §[}]). On the other 
hand, for M > 10 _1 - 7 Mq T c is below the CMB temperature 
and the solid phase warms the cloud rather than cooling it. 

3. THERMAL EQUILIBRIUM 

In thermal equilibrium, the cloud temperature is set by the 
balance between cooling and heating. We assume that clouds 
in the Galactic halo are heated primarily by cosmic rays. The 
local interstellar cosmic -ray ionisation rate in the Galactic disc, 
- 3 x lO-^s" 1 ^ 1 (Webber 1998), implies a heating rate 
T - 3 x 10" 4 ergg- 1 s" 1 (Cravens & Dalgarno 1978). The 
cosmic-ray heating in the halo is uncertain but should be some- 
what lower, say ~ 10~ 5 ergg _1 s _1 . In Fig. 0, we show the 
cooling rates for clouds of masses 10 -7 — 10 _1 MVIq. The solid 
curves are contours of constant optical depth; the dashed curve 
shows the optically-thick limit and represents the maximum 
cooling rate for each cloud mass. It appears that solid hydro- 
gen can provide the necessary cooling for planetary-mass gas 
clouds at the cosmic -ray heating rates expected in the Galactic 
disk and halo. 

The upper panel of Fig. [TJ shows that there are typically three 
equilibrium temperatures available for cloud masses between 
10~ 7 5 and 10 _17 Mq at the expected heating rates: solid H 2 
provides one barely above T b and one a few degrees higher; the 
gas-phase coolants provide an equilibrium above 30 K. We now 
show that thermal stability requires A to be a decreasing func- 
tion of T, and therefore only the second of these three equilibria 
is stable. 

In virial equilibrium the total energy per unit mass is approx- 
imately — |fcT/ [i (the internal excitation of the gas is negligible 
at the low temperatures of interest here), so the thermal evolu- 
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FIG. 2. — Solid H2 cooling rate versus cloud mass for different optical depths. The curves, labelled by t, are truncated at high masses when the solutions become 
thermally unstable. The optically thick limit, corresponding to black-body emission, is labelled "00". The dashed curve delineates where half of the hydrogen is in 
solid form, an approximate upper limit for thermal stability. 



tion of the cloud is determined by 



3k dT 
2fl~dJ 



A-r. 



(8) 



In the absence of heating the cloud contracts on the Kelvin- 
Helmholtz time-scale ^kh = §£;T/(/iA). Note that this time- 
scale can be a substantial fraction of the Hubble time for tem- 
peratures of a few Kelvin: the dashed curves in figure |l] show 
the cooling rate that yields £kh = 10 Gyr. In thermal equi- 
librium cosmic -ray heating replaces the energy radiated away 
by the cloud implying, for example, £kh ~ 2 x 10 6 yr for 
r ~ 1CP 5 erg g _1 s _1 (at T ~ 10 K). This is much greater 
than the sound crossing time (~ 10 2 yr), so the response of 
a cloud to dynamical perturbations is adiabatic, to a good ap- 
proximation, and dynamical stability is assured. However, eq. 
(^|) shows that perturbations to the cloud temperature grow or 
decay as e at where 



« = t KH?4(A-r) 



A dT' 



(9) 



and the right-hand side of this equation is evaluated at the equi- 
librium temperature. Thus a cloud is thermally stable only if a 
decrease (increase) in cloud temperature leads to cooling out- 
stripping (lagging) heating. For cosmic ray heating, T is inde- 
pendent of T if the column through the cloud is insufficient to 
cause significant attenuation of cosmic rays (changes in tem- 
perature affect the cloud's column density through the virial re- 
lationship R oc 1/T). Thermal stability then requires that A 
be a decreasing function of T, and we conclude that only the 
equilibrium on the high-temperature shoulder of the solid hy- 
drogen cooling curve is stable. In fact the column density of 
each cloud (~ 10 2 gem -2 : Walker 1999) is sufficient to stop 
sub-GeV cosmic -ray protons (and all electrons), leading to a 
dependence of T on T; this dependence is too weak to affect 
our conclusions concerning stability. 

4. DISCUSSION 

The suggestion that cold gas could comprise a significant 
fraction of the Galaxy's dark matter is not new, previous pro- 
posals include: a fractal medium in the the outer reaches of 
the Galactic disk (Pfenniger et al. 1994); isolated halo clouds 
(Gerhard & Silk 1996); and mini clusters of clouds in the halo 



(de Paolis et al. 1995; Gerhard & Silk 1996). However, to date 
there has been no compelling reason to believe that isolated, 
cold gas clouds - as inferred by Walker & Wardle (1998) - 
could support themselves for long periods against gravitational 
collapse. We have shown that such clouds can be stabilised by 
the precipitation/sublimation of particles of solid hydrogen (or 
by the condensation/evaporation of droplets of liquid hydrogen) 
if these particles dominate the radiative cooling of the cloud. 
The key feature which confers thermal stability is that these 
particles are destroyed, hence cooling becomes less efficient, as 
the cloud temperature increases. This feature will be present in 
any model where condensed hydrogen is the principal coolant, 
and consequently we expect that more sophisticated structural 
treatments will also admit stable solutions. 

The masses of thermally stable clouds lie in the approxi- 
mate range 1O~ 7 ' 5 -1O~ 1 ' 7 M0. The lower limit is increased to 
10 _6 Mq if subject to cosmic ray heating similar to that in the 
Galactic disc for an interval J> kT/T ~ 10 5 yr. As halo clouds 
take much longer than this to pass through the cosmic -ray disc, 
this limit is appropriate even for a halo cloud population. For 
cloud masses in the range 10~ 6 -10~ 17 Mq the radiative cool- 
ing simply readjusts, on the timescale £kh> to maintain equilib- 
rium as r varies through the orbit (~ 10 8 yr) of a cloud around 
the Galaxy. 

The typical particle radius a is constrained by the require- 
ment that the clouds not produce significant extinction at opti- 
cal wavelengths: the geometrical optical depth of the particles, 
T g w T/C m ap s should be less than 1. Adopting T = 5K 
and t — 0.01 this translates to a(mm) J> (10~ 4 cm)/A2. 
Millimeter-size particles settle to the centre of the cloud in 
~ 10 4 yr, this time is shortened if A2 is significantly less than 
10~ 4 cm and the particles are required to be larger. Settling 
may be counteracted by convective motions or sublimation re- 
sulting from the higher temperatures deeper within the cloud 
— issues that must await a more sophisticated treatment of the 
cloud structure. 

If the hypothesised population of clouds exists, their ther- 
mal microwave emission may be detectable as a Galactic con- 
tinuum background at temperatures just above the Cosmic Mi- 
crowave Background; a Galactic component of this kind has in 
fact been isolated in the COBE FIRAS data (Reach et al. 1985). 
One would like to compare these data with the theory presented 
here, but it is difficult to predict the total microwave intensity 
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for our model because the distribution of cosmic-ray density 
away from the Galactic plane is only loosely constrained (see 
Webber et al. 1994). A similar uncertainty afflicts the modelling 
of 7-ray production from baryonic material in the Galactic halo 
(cf. de Paolis et al. 1995; Salati et al. 1996; de Paolis et al. 
1999; Kalberla, Shchekinov & Dettmar 1999). Nevertheless, 
microwave and 7-ray emissivities are each proportional to the 
local cosmic -ray flux (assuming the cosmic -ray spectrum does 
not vary greatly), so we can write I M ~ I^j^l ' j-y, for emissiv- 
ities j and intensities /. At high latititudes the Galactic 7-ray 
background is 1 1 ~ phcm -2 s _1 sr~\ above 1 GeV 

(Dixon et al. 1998); local to the Sun the corresponding (opti- 
cally thin) emissivity is 1.1 x 10~ 3 ph s _1 g _1 sr _1 (Bertsch 
et al. 1993). Thus for a cosmic-ray heating rate (again, local to 
the Sun) of T = Airj^ - 3 x 1CT 4 ergs" 1 g" 1 (Cravens & Dal- 
garno 1978), we expect 1^ ~ 2xl0~ 8 erg cm -2 s _1 sr _1 . This 
is roughly 1 % of that observed in the FIRAS cold component 
at high latitude (Reach et al. 1985), so the microwave data do 
not exclude the possibility of a cold cloud population heated by 
cosmic rays. Sciama (1999) proposed that all of the cold excess 
may be accounted for by cosmic-ray heating of cold clouds, but 



this appears to be based on an overestimate of the gamma-ray 
flux, and an underestimate of the high-latitude FIRAS flux. 

5. CONCLUSIONS 

We have demonstrated that, by virtue of the solid/gas 
phase transition of hydrogen, cold, planetary-mass Galactic gas 
clouds can be thermally stable even when they are heated at 
a temperature-independent rate. Our analysis applies to the 
present epoch, with the microwave background temperature at 
Tb < 3 K; for background temperatures ;> 6 K, our model 
admits no stable mass range. Consequently the longevity of the 
clouds at redshifts z J> 1 is problematic. We cannot, however, 
hope to address this issue until a firm theoretical basis for the 
formation of such clouds has been established. 
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